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The rupture of dry frictional interfaces occurs through the propagation of fronts breaking the 
contacts at the interface. Recent experiments have shown that the velocities of these rupture fronts 
range from quasi-static velocities proportional to the external loading rate to velocities larger than 
the shear wave speed. The way system parameters influence front speed is still poorly understood. 

Here we study steady-state rupture propagation in a one-dimensional (ID) spring-block model of 
an extended frictional interface, for various friction laws. With the classical Amontons-Coulomb 
friction law, we derive a closed-form expression for the steady-state rupture velocity as a function 
of the interfacial shear stress just prior to rupture. We then consider an additional shear stiffness 
of the interface and show that the softer the interface, the slower the rupture fronts. We provide 
an approximate closed form expression for this effect. We finally show that adding a bulk viscosity 
on the relative motion of blocks accelerates steady-state rupture fronts and we give an approximate 
expression for this effect. We demonstrate that the ID results are qualitatively valid in 2D. Our 
results provide insights into the qualitative role of various key parameters of a frictional interface on 
its rupture dynamics. They will be useful to better understand the many systems in which spring- 
block models have proved adequate, from friction to granular matter and earthquake dynamics. 


I. INTRODUCTION 

Extended frictional interfaces under increasing shear 
stress eventually break and enter a macroscopic sliding 
regime. They do so through the propagation of a rupture¬ 
like micro-slip front across the whole interface. The prop¬ 
agation speed of such fronts is typically of the order of 
the sound speed in the contacting materials, which made 
them elusive to measurements until the rise of fast cam¬ 
era monitoring of frictional interfaces in the late 1990s. 
It is now well established that a whole continuum of front 
propagation speeds v c can be observed along macroscopic 
frictional interfaces, from intersonic (between the shear 
and compression wave speeds of the materials, c s and 
Ck, see e.g. [1]) to quasi-static (proportional to the ex¬ 
ternal loading rate, see e.g. [2, 3]), going through sub- 
Rayleigh fronts (y c < cr with cr the Rayleigh wave 
speed, see e.g. [4-6]) and slow but still dynamic fronts 
(v c ~ 0.01 to O.Icr, see e.g. [4]). This huge variety in 
observed speed triggered the natural question of what 
the physical mechanisms underlying speed selection of 
micro-slip fronts are. 

Experimentally, it has been shown that the larger the 
local shear to normal stress ratio r/p just prior to rupture 
nucleation, the faster the local front speed [5]. This re¬ 
sult is consistent with the observation that a larger shear 
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stress promotes intersonic rather than sub-Rayleigh prop¬ 
agation [7, 8]. The observed relationship between pre¬ 
stress and front speed has been reproduced in simulations 
for both fast [9-11] and slow [11] fronts. On other aspects 
of front speed, models are ahead of experiments, with var¬ 
ious predictions still awaiting experimental verification. 
Among these are the following: (i) Models with simple 
Amontons-Coulomb (AC) friction [9, 11-13] have shown 
that front propagation speed is controlled by f = > 

with n s and p ,k the local static and kinematic friction 
coefficients, t thus appearing as a generalisation of the 
parameter r/p used to analyse the experimental data, 
(ii) A model with velocity-weakening AC friction [10] has 
suggested that front speed is direction-dependent, with 
different speeds for fronts propagating with and against 
the shear loading direction, (iii) Two-dimensional (2D) 
spring-block models [11, 13] and ID continuous models 
[14, 15] have shown that front speed Vf ron t is proportional 
to some relevant slip speed w s i; p , with a relationship of 

j-u j__ ry-vshear modulus of the interface,,, 
une type Ufront ~ stress drop during rupture ^ sli P* 

Giving quantitative predictions of front speed is diffi¬ 
cult for at least two reasons. First, any real interface is 
heterogeneous at the mesoscopic scales at which stresses 
can be defined (scale including enough micro-contacts), 
due both to intrinsic heterogeneities of the surfaces and 
to heterogeneous loading. Thus, even if the front speed 
was selected only locally, i.e. as a function of the lo¬ 
cal stresses and local static friction threshold, the front 
speed would still be varying with front position along the 
interface. Second, models actually show that front speed 
can have long transients [13] (extending over sizes com- 
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parable to that of the samples used in a number of exper¬ 
iments), even for carefully prepared homogeneous inter¬ 
faces. This suggests that the instantaneous front speed 
does not only depend on local quantities, but rather on 
the slip dynamics along all the broken part of the inter¬ 
face and from all times after front nucleation. Here, we 
overcome these difficulties by (i) considering fronts prop¬ 
agating along homogeneous interfaces and (ii) focusing 
on front speed only in steady state propagation, i.e. af¬ 
ter transients are finished. For the sake of simplicity and 
to enable analytical treatment, we use a ID spring block 
model for the shear rupture of extended frictional inter¬ 
faces, introduced in [ 12 ] to study the propagation length 
of precursors to sliding [9, 16-20]. Whereas the dynam¬ 
ics of multiple successive events in the macroscopic stick- 
slip regime of this model was discussed in detail in [ 12 ], 
here we focus on the steady state propagation of a single 
rupture front. We show in the discussion that the main 
results obtained using the ID model still hold in a 2D 
extension of the model. 

We emphasise that fully dynamic (as opposed to cellu¬ 
lar automata) spring-block (or spring-mass) models have 
previously been widely used in the literature to model 
not only friction [see e.g. 9, 11, 12, 17, 21, 22] and earth¬ 
quake dynamics [see e.g. 23-25], but also, among others, 
self-organized criticality in nonequilibrium systems with 
many degrees of freedom [e.g. 26], adsorbed chains at sur¬ 
faces [e.g. 27], fluctuations in dissipative systems [e.g. 28] 
or creep in granular materials [29]. 

Rupture velocities in ID spring-block models have 
been studied previously [30-32] in the framework of the 
Burridge-Knopoff (BK) model [23]. In the BK model, a 
chain of blocks and springs is loaded uniformly from the 
top through an array of springs connected to a rigid rod. 
Note that this loading configuration differs from the one 
used in the present paper, in which the chain of blocks is 
loaded from one edge. In [30], the rupture speed of the 
BK model with velocity weakening friction was obtained 
in the case of a uniform loading exactly at the local slip¬ 
ping threshold. The model was found to have a range 
of possible propagation velocities among which one is se¬ 
lected dynamically. The rupture velocity was also found 
to be resolution-dependent. This resolution problem was 
solved in [31] by introducing a short-wavelength cutoff, 
obtained by adding Kelvin viscosity to the model. Rup¬ 
ture velocities in the BK model with Amontons-Coulomb 
friction were studied in Muratov [32], and found to have 
a unique solution for any given value of the initial shear 
stress at the interface, with a well-defined continuum 
limit. We compare our results to those of [30-32] in Sec¬ 
tion IV. 

The paper is organised as follows: We first describe our 
model and derive its non-dimensional form (Section II). 
We then present our results for the velocity of steady- 
state front propagation as a function of the pre-stress on 
the interface prior to rupture (Section III), for three vari¬ 
ants of the model. We start with a simple AC friction law 
and obtain a closed form equation for the front velocity. 


We then add either a bulk viscosity or an interfacial stiff¬ 
ness, and provide for each an approximate equation for 
front speed. In Section IV, we discuss our results in the 
light of a 2D model. Conclusions are in Section V. Four 
appendices provide additional mathematical details. 


II. MODEL DESCRIPTION 

We investigate the propagation of micro-slip fronts in 
the ID spring-block model originally introduced by Mae- 
gawa et al. [17] to study the length of precursors to slid¬ 
ing. It has been later improved by us [12] to include 
a bulk viscosity and a friction law allowing for a finite 
stiffness of the interface. A schematic of this minimalis- 
tic model is given in Fig. 1. The slider is modelled as a 
chain of blocks with mass m = M/N connected in series 
by springs with stiffness k = (N — 1 )ES/L, where M is 
the total mass of the slider, N is the number of blocks, E 
is the Young’s modulus, S is the cross-section area and 
L is the length of the slider. The applied normal force 
on each block n is given by p n = F^/N, where Fn is 
the total (uniformly) applied normal force. The tangen¬ 
tial force Ft is applied at the trailing edge of the system 
through a loading spring with stiffness K. One end of 
this spring is attached to the trailing edge block (block 
1 ), while the other end moves at a (small) constant ve¬ 
locity V. 



FIG. 1: Schematic of the model system. The slider is 
modelled as N blocks with mass m connected by 
springs of stiffness k. The trailing edge block (block 1) 
is slowly driven through a loading spring of stiffness K. 
Each block n is also submitted to a normal force p n , a 
friction force f n and a viscous damping force F%, which 
is described by the viscous coefficient 77 


The equations of motion are given by 

mu n = F% + F% + f n , 1 < n < N, (1) 

where u n = u n ( t ) is the position of block n as a function 
of time relative to its equilibrium position (in the absence 
of any friction force) and " denotes the double derivative 
with respect to time t. The forces F 1 *, F r r j and f n are the 
total spring force, relative viscous force and the friction 
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force on block n, respectively, and are given by 

J k{u 2 - Mi) + Ft, n= 1 

k(u n+1 -2u n + u n - 1 ), 2 < n < N — 1 (2) 

k(u N _ 1 -u N ), n = N, 

{ t](u 2 -ui), n = 1 

T](u n+ 1 - 2u n + u n -i), 2 < n < -/V — 1 (3) 

ri(ii N -i - un), n = N, 

with the tangential load Ft given by 

Ft = K (Vt — Mi). (4) 


Next, we introduce our dimensionless variables, u n = 
u' n /U for block positions, t = t/T for time, and x = x/X 
for horizontal positions. Substituting these back into 
Eqs. (1) and ( 6 ) to ( 8 ), yields our dimensionless equa¬ 
tions of motion. We make the following choices for the 
scaling parameters U, T and X: 

U= (Ms T = yfrnjk, X = a, (10) 

k 

where and ^ are the static and kinetic friction coef¬ 
ficients in AC-like friction laws and a = L/{N — 1). The 
dimensionless equations of motion become 


In the following we will simply use the term “viscosity” 
to mean the viscous coefficient 77 . We consider two dif¬ 
ferent functional forms for the friction force f n , one cor¬ 
responding to the rigid-plastic-like Amontons-Coulomb 
(AC) friction law, discussed in Section III A, and one to 
the elasto-plastic like friction law introduced in Amund¬ 
sen et al. [ 12 ] allowing for a finite stiffness of the interface, 
discussed in Section IIIB. 

Before solving Eq. (1) it is instructive to rewrite it on a 
dimensionless form to derive the combination of param¬ 
eters controlling the behaviour of the system. Here we 
derive the dimensionless equations of motion for a generic 
friction force f n and will later consider the two special 
cases discussed above, (i) AC friction in Section III A 
and (ii) with tangential stiffness of the interface in Sec¬ 
tion III B. 

We begin by eliminating the initial positions of all 
blocks, u n (0), from the block positions u n (t). Any move¬ 
ment can be described by u' n {t) defined by 

u n (t) = u n ( 0) + u' n (t), (5a) 

Unit) =u' n (t), (5b) 

Unit) = Unit), (5c) 


u n — F n + F% + fn, (11) 


with 


Fn 


pv 


fn 


Ft 


u 2 — u\+ Ft, n = 1 
U n+ 1 — 2 U n + U n -!, 2 < 71 < N — 1 
un-i — un, n = N, 

! jj(it 2 - fti),. _ n= 1 

tjiUn+l - 2u n + Un- 1), 2 < U < N - 1 
f)iu N -!-u N ), n = N, 

1 If , - \ _ T~n F fn 

777 \Jn ~T T n ) — ~ \ 5 

kU (/i s fJ'kjPn 


—~ K {vt - 7) = i<(vi- no, 

mu 


( 12 ) 

(13) 

(14) 

(15) 


where ' now denotes the derivative with respect to i and 
not t. Note that for convenience, and being of frictional 
origin, the initial shear force r„ has been included in 
the effective dimensionless friction force f n in Eq. (14) 
rather than in F n ■ The dimensionless relative viscosity 
is defined as 


V = 



(16) 


i.e. the position of a block is the position it had at t = 0 
plus any additional movement u' n (t). The forces F%, Ff> 
and Ft then become 


. [ 

f H u 2 - 

7) + 

T\ + F^, 

n = 

1 


II 

fef 

, Hu'n+l 

-2 u'n 

+ u'n- 1 ) + T n , 

2 < 

n < N — 

1 

1 

{ Hu'n _1 

- u' N 

) + tn, 

n = 

N, 








(6) 

I 

[ viu'i - 

7), 

n = 

1 



II 

viu'n+l 

-2 U' n 

+ u'n-i), 2 < 

VI 

e 

N — 1 

(7) 

I 

[ viu ' N -1 

~u’ N 

), n = 

N, 



= K(Vt - u‘ 

i), 




(8) 


and we have introduced 


V = 


V 

ujf’ 


(17) 


The velocity of sound in this model is given by [33] 


i> s = a\l —, 
m 


(18) 


and in our dimensionless units this becomes 


Sa= = 


(19) 


where we have introduced a new force r„, the initial shear 
force, given by 

f k{u 2 (0) — Ui(0)) — Kui{0), n= 1 
r n = < k(u n+ i(0) - 2u n (0) + u„_!(0)), 2 < n < — 1 

[ A:(ujv-i(0) - u N (0)), n = N. 

( 9 ) 


which was the reason for our choice of X. 

Looking at our new dimensionless set of equations it is 
clear that the number of parameters has been reduced. In 
addition to the dimensionless friction force f n , only K, V 
and fj, the dimensionless driving spring constant, driving 
velocity and relative viscosity, respectively, will impact 
the evolution of the dimensionless block positions. 
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III. STEADY STATE FRONT PROPAGATION 

Most previous studies of rupture front propagation in 
spring-block models have initialised the models with the 
shear stresses set to zero. External loading is then ap¬ 
plied, and the evolution of the systems in time is stud¬ 
ied [9, 12, 17- 19, 21]. The interface states at the time 
of nucleation of micro-slip fronts are selected naturally 
through the evolution of the system in time, mimicking 
the experimental setups of e.g. [16, 17]. 

To facilitate systematic study of the front velocity in 
our system, as previously done in e.g. [11, 13], we prepare 
a desired interface state at the time of rupture nucleation 
and look at the resulting rupture dynamics. The initial 
state is governed by the initial shear stresses, t„, which 
is one of the important parameters in the effective fric¬ 
tion force f n . In this paper we discuss steady-state front 
propagation, i.e. fronts propagating at a constant veloc¬ 
ity v c , and for this reason all surface properties, including 
T n , are kept homogeneous throughout the interface. For 
convenience, the block index n will therefore be dropped 
in the following discussion where possible. We also make 
the assumption that the driving spring constant K is 
much smaller than the material spring constant k. i.e. 
K -C 1, which means Ft can be treated as constant dur¬ 
ing the front propagation. This assumption is valid both 
in models studied previously (e.g. [9, 12, 17]) and in ex¬ 
perimental studies (e.g. [5, 17]). 

Here we use both numerical and analytical tools to 
study steady-state rupture fronts in ID spring-block 
models. In Section III A we measure the front speed in 
our model with AC friction. We compare the cases with¬ 
out and with bulk viscosity fj and show that in a few 
special cases closed-form expressions for the front veloc¬ 
ity as a function of model parameters may be obtained. 
In Section III B we study the impact of a finite stiffness 
of the interface before concluding with some remarks on 
the complete model. 


A. Amontons—Coulomb friction 


Perhaps the simplest dry friction law in wide-spread 
use is Amontons-Coulonrb friction, which introduces 
static and dynamic friction coefficients /x s and /.tk, re¬ 
spectively. We impose this law locally on each block in 
our system as in [9, 12, 17], i.e. a block has to overcome 
a friction threshold p s p to start sliding, during which it 
experiences a force p^p. The friction force /„ is therefore 
given by 


f < Ms p, Un = 0 

\ —Sgn(u n )/XkP, Un ^ 0 , 


( 20 ) 


where, when u n = 0, f n balances all other forces acting 
on block n. Blocks are assumed to repin to the track 
when their velocity becomes 0 and will only start moving 
again if the sum of all forces, except the friction force, 
again reaches the static friction threshold p s p. 


For sliding blocks we insert Eq. (20) into Ecp (14) and 
obtain the effective dimensionless friction force f n : 


j± = t/p T Mk 
Ms - Mk 


( 21 ) 


Note the necessary separation into f + and / _ at this 
point, where / + applies if the block is moving in the 
positive direction and f~ if it is moving in the negative 
direction. This is related to the change in sign of the 
friction force f n as the block velocity changes between 
being positive and negative. 

In the model, a front propagates in the following way: 
The driving force increases on block 1 up to the local 
static friction level. As block 1 moves, the tangential 
force on block 2 increases, eventually reaching its static 
friction threshold, and starts to move. We interpret the 
successive onset of motion of blocks as the model equiva¬ 
lent of the micro-slip fronts observed in experiments and 
define the local front velocity as the distance between two 
blocks divided by the time interval between the rupture 
of two neighbouring blocks. We label these two blocks 
n = i and n = i + 1, and denote the time between the 
onset of motion of these two blocks by At,. Since the 
material springs are very stiff, the distance between two 
neighbouring blocks can be approximated to be a, inde¬ 
pendent of time. We define the rupture velocity v c as 


v c 


a 

At? 


V c 



Vc 

7 

V s 


( 22 ) 


where we have used the velocity of sound in Eq. (18). 

A block begins to move if the forces on it reach the 
static friction threshold, 

F* + FZ = ±p s p, (23) 

which in dimensionless units becomes 


F n + F^ = ±l-f ± . (24) 

A rupture initiates when the total force on block 
1 exceeds the static friction threshold, while all other 
blocks are still unaffected. We initialise the system with 
t + F! t = m s P, i-e. 

f + Ft = 1, (25) 


where 

- = t /p - Mk = j+ 

Ms - Mk 


(26) 


Note the definition of f, which we will later show to be 
a very important model parameter. We exclusively con¬ 
sider positive initial shear forces, the maximum being 
restricted by the static friction threshold. Consequently, 
all values of t, the dimensionless initial shear force, lie 
between — Mk/(Ms — Mk) and 1- 

To summarise, the equations of motion for the system 
are given by Eqs. (11) to (15) and (21), and rupture 
initiates when Eq. (25) is satisfied. We will now proceed 
by first considering the simplest case where fj = 0, before 
studying the effect of introducing a bulk viscosity. 
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1. Model without bulk viscosity (r) = 0) 


As in the model by Maegawa et al. [17], we first take 
fj = 0 and apply AC friction locally at each block. Since, 
as discussed above, we keep Ft constant during rupture, 
and since its initial value is given by the rupture criterion 
(Eq. (25)) only two parameters remain that control the 
front propagation, f = j + and /“. We first want to 
identify for which values of steady-state ruptures can 
be supported. 

We have tried many different values for r, some exam¬ 
ples of which are shown in Fig. 2a. We observed that 
blocks never move in the negative direction for as long 
as the front propagates, which means that f~ becomes 
irrelevant, and we are left with only one parameter, f, 
controlling front propagation. We have found that, for 
a steady-state rupture to be supported, f > 0 is re¬ 
quired. The natural restriction r/p < p s places another 
constraint on f, and we can conclude that steady-state 
ruptures occur only if 

0 < f < 1. (27) 


The straightforward way to compute the steady-state 
rupture velocity as a function of f is to solve Eqs. (11) 
to (15) and (21) explicitly in time for a given value of f. 
Fronts go through a transient, as seen in Fig. 2a, before 
reaching the steady-state velocity. As shown in Fig. 2b, 
the transient length is strongly dependent on the pre¬ 
stress, f, and extrapolation is necessary to estimate the 
final steady-state rupture velocity. We extrapolate the 
rupture velocity by fitting a first order polynomial to the 
curve v c (l/n) for the last few (~ 50) blocks towards the 
leading edge. These extrapolated steady-state velocities 
are plotted as a function of t in Fig. 3 as red circles. 

Alternatively, equations for the steady-state front ve¬ 
locity can be derived from Eqs. (11) to (15) and (21). We 
provide this derivation in Appendix A 1, with Eqs. (A3), 
(A5), (A6) and (A8) the final set of equations. The nu¬ 
merical scheme used to solve these equations is detailed 
in Appendix B 2. We show in Fig. 3 the steady-state front 
velocity v c as a function of f obtained by solving these 
equations numerically as blue crosses. This solution is 
seen to match the extrapolated front velocities obtained 
previously. 

It is also possible to solve the equations for the steady- 
state rupture velocity, Eqs. (A3), (A5), (A6) and (A8), 
analytically using an iterative approach. The solution 
technique is identical to the one used by Muratov [32], 
and we present the detailed calculation in Appendix B 1. 
The result is a series expansion of f as a function of 
z = l/va given in Eq. (B8). In the present case with 
fj = 0, we get 


_ _ z 2 z 4 z b 5z® 

T ~ 2 8 16 _ 128 


10 


7z 

256 


+ O (z 12 ) , (28) 


which we recognise as the series expansion of 

-+0 (z 12 ) . (29) 


n -j z 2 z 4 z 6 5z 8 7z 10 

^ ~~2 8~ ~ 16 ~ 128 _ ~256 



Normalised system length 


(a) Simulated front velocity as a function of position along the 
interface for several different values of f = obtained with 
N = 200. From bottom to top, f = —0.1, —0.01, 0.01, 0.1, 0.5. 
Velocities are seen to approach a steady-state velocity for f > 0, 
while for f < 0 the fronts slow down and eventually arrest. 



T 

(b) Transient length, defined as the block number where the front 
velocity reaches 97 % of the steady-state value, normalised by 
system size as a function of f obtained with N = 1000. 

FIG. 2: Colour online. Rupture velocity as a function of 
position along the interface (a) and front velocity 
transient lengths (b) for various initial shear stress. 


Consequently, we have a closed-form expression for the 
front velocity as a function of the initial shear stress: 

f = y/l- Vc 2 or v c = 1 . (30) 

V 1 — T 2 

This solution is plotted as the black solid line in Fig. 3, 
which matches the numerical results perfectly. 

To summarise, the steady-state front velocity v c in 
the model with AC friction and fj = 0 only depends 
on the dimensionless initial shear stress, r, and is given 
by Eq. (30), plotted in Fig. 3. The front velocity in¬ 
creases with increasing f , as expected. All steady-state 
front velocities are supersonic, with values approaching 
the sound velocity as f —► 0 and infinity as t —i 1. The 
latter is easily explained: as f —i 1, every block will be 
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FIG. 3: Colour online. Steady-state front velocity v c as 
a function of initial stress, t. v c approaches 1 as f —> 0 
and infinity as f —i 1. Red circles are the extrapolated 
steady-state front velocities from simulations (as in 
Fig. 2a), blue crosses are the numerical solution of the 
steady-state equations (Appendix Al), and the solid 
black line is the analytical solution given by Eq. (30). 



T 


FIG. 4: Colour online. Steady-state front velocity v c as 
a function of initial shear stress, f, for several different 
values of the relative viscosity fj. From bottom to top, 
fj = 0 (blue, circles), \/07T (green, crosses), 1 (red, 
pluses) and \/l0 (cyan, squares). The point markers are 
the numerical solutions of the steady-state equations, 
and the solid lines are the semi-empirical closed-form 
expression in Eq. (36). 


infinitely close to the static friction threshold, and in¬ 
finitesimally small movement of the neighbouring blocks 
is enough to set them into motion. The time between 
the triggering of neighbouring blocks will therefore ap¬ 
proach 0, causing the front velocity to approach infinity, 
see Eq. (22). This is a known feature of spring-block 
models [34], We discuss these results in more detail in 
Section IV. 

This model and the above results serve as our reference 
for investigating now the effect of a bulk viscosity and a 
tangential stiffness of the interface on the steady-state 
front velocity. 


2. The effect of a bulk viscosity 

In this section we study the effect of the relative viscos¬ 
ity fj , identical to the one used in Amundsen et al. [12] 
to smooth grid-scale oscillations during front propaga¬ 
tion [31, 35]. Physically it is a simple way of introducing 
energy dissipation that will occur during deformation. As 
in Section III A 1 we have found that steady-state rup¬ 
tures may occur if 


0 < f < 1, (31) 

independent of the value of the viscosity, fj. Similarly we 
have also found that blocks exclusively move in the posi¬ 
tive direction as the front propagates along the interface. 
Consequently, we have two parameters controlling front 
velocity in this system, f and fj. 

The steady-state equations are solved numerically as 
in Section III A 1 for several different values of fj, and 
we plot v c (t) in Fig. 4. As in the model with fj = 0, the 
steady-state equations can be solved analytically, and the 


result is a series expansion of f in terms of z = 1 /v c and 
fj. For brevity we do not reproduce this expansion here, 
it is given in Eq. (B8). 

Unfortunately we have been unable to find a closed- 
form expression for the front velocity as a function of f 
and fj from Eq. (B8). The special case fj = 1, however, 
yields 

f = 1 - z + O (z 12 ) , (32) 

i.e. 

f = l—— or v c = —. (33) 

v c 1 -t 

We use this to derive a semi-empirical expression for 
t(v c ). From Eqs. (30) and (33), we can write 

/ 1 — t; -1 \ 1/2 

t{v c , rj = 1) = t(v c , rj = 0) _ c _ 1 ) . (34) 

V1 + v c J 


The factor [(1 — ^)/(l -F ^r)] 1 ^ 2 can be thought of as a 
more general scaling factor including the fj dependence, 
but taken at fj = 1. We have found that we can estimate 
the fj dependence with an exponent of fj/ 2, which fits the 
numerical solution fairly well. We therefore propose the 
following approximation 


v(v c ,il ) 


t(v c ,t} = 0) 


f 1 - w 1 

\ 1 T v c 


fj/2 


Vl - Vc 2 


( 1-Vc 1 
\ 1 + Vc 1 


fj/2 


(35) 

(36) 


We plot Eq. (36) in Fig. 4, and it is seen to match the 
numerical solution perfectly for fj = 0 and fj = 1 as ex¬ 
pected, and to yield acceptable accuracy for 0 < fj < 1. 









7 


The value of fj = \/0T adopted by Amundsen et al. [12] 
is within this range. The accuracy of the semi-empirical 
expression deteriorates for fj > 1, which corresponds to a 
regime where all waves are overdamped [12]. 

In Fig. 4 the relative viscosity is seen to not alter the 
limiting behaviour as f —> 0 and f — > 1, for which the 
front velocity still approaches v c = 1 and v c —> oo, respec¬ 
tively. An overall increase of front velocities compared to 
the fj = 0 case, discussed in Section III A 1, is, however, 
seen. 

The reason for this particular behaviour is that the 
relative viscosity serves to dampen out (reduce) relative 
motion between blocks. At the front tip, the rightmost 
moving block is increasing the load on the leftmost stuck 
block through the material spring connecting them and 
through the relative viscous force. This viscous force acts 
in the direction of movement of the moving block. This 
causes the static friction threshold to be reached sooner, 
and the stuck block starts moving earlier than it would 
have with a smaller value of fj. Note that the stuck block 
will act with an equal and opposite force on the moving 
block, slowing it down. This effect remains small, how¬ 
ever, due to the large momentum of the blocks behind the 
front tip. Overall, Ai n . the time between the rupture of 
block n and n + 1, is reduced. 

B. Elasto-plastic like friction law 

As discussed in Amundsen et al. [12], in the model 
with AC friction, only the first block will experience the 
tangential loading force. This causes an unphysical res¬ 
olution dependence in the model, which was improved 
considerably by introducing a finite tangential stiffness 
of the interface. In addition, the interface between the 
slider and the base is indeed elastic (see e.g. [2]), a feature 
which is often accounted for in models using an ensemble 
of interface springs to model the micro-contacts binding 
the slider and base together [11, 13, 21, 36, 37]. 

We introduce a tangential stiffness of the interface as in 
Amundsen et al. [12] by modifying the static Amontons- 
Coulomb friction law to include springs between the 
blocks and the track. Each block bears one interface 
spring having a breaking strength equal to the static 
friction threshold p s p. When a spring breaks, dynamic 
friction p^p applies until the spring reattaches when the 
block velocity becomes zero. The spring is reattached 
such that at the time of reattachment the total force on 
the block is zero. 

In this section we study front velocity as a function 
of f and the stiffness of the interface springs, k t . For 
attached blocks the friction force is given by 

fn = -h ( U n (t) - <‘ ick (t)) 

= ~kt ( u' n (t) - < ick '(f)) - h (u n ( 0) - < tick (0)), 

(37) 

where w^ lck (t) is the position of the attachment point of 


the spring. At t = 0 the total force on all blocks is zero, 
i.e. f n (t = 0) = T n , and we get 

fn = -kt (u' n {t ) - < tIck '(t)) - T n . (38) 

Using Eq. (14) we obtain an expression for the dimen¬ 
sionless friction force, 



Ms Pn Mk Pn 


-k t (<(f)-< tick '(f)) 

=- ^(40) 

Ms Pn Mk Pn 

= -~k (unit) - < tick (t)), (41) 

where k = kt/k. 

The rupture criterion is modified as it is a condition 
on the strength of the interface springs. It is given by 

kt (u n {t) - < tlck (f)) = p s p, (42) 

which in dimensionless variables becomes 

k ( Unit ) ~ A^ lck (i)) + Ai = 1. (43) 

As discussed above, interface springs reconnect when the 
block velocity becomes zero and reconnect at zero total 
force: 

0 = F n + K + f n (44) 

= F n + F* - k (■ u n (t ) - < tick (f)) , (45) 

which yields 

p /y.stick\ _i_ pfj /y.stick\ 

-stick Q* tick) = u n (t stick ) - -Alb - _ - 2, (46) 

k 

where the only new parameter introduced is the dimen¬ 
sionless interface stiffness k = kt/k, and £ stlck is the time 
at which the block velocity becomes zero. w^ tlck (f) stays 
constant for t > t stlck until the block reattaches after 
another detachment event. 

For simplicity we investigate the behaviour of this 
model without the relative viscosity here, fj = 0, but 
this assumption is relaxed in Section IIIC. As in Sec¬ 
tion III A 1 we have investigated when steady-state rup¬ 
tures occur and found that it is again restricted to 

0 < r < 1, (47) 

independent of the value of the interface stiffness, k. Sim¬ 
ilarly we have also found that blocks exclusively move in 
the positive direction as the front propagates along the 
interface. Consequently, we have two parameters con¬ 
trolling the rupture velocity in this system, t and k. 

We solve the steady-state equations (derived in Ap¬ 
pendix A 2) and show in Fig. 5 the front velocity as a 
function of the dimensionless initial shear force t for var¬ 
ious values of the interface stiffness k. In the limit k —> oo 
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FIG. 5: Colour online. Steady-state front velocity v c as 
a function of initial shear stress, t, for several different 
values of the interface stiffness k. From top to bottom, 
k = 10 3 , 10 1 , 1, 10 _1 , 10 -3 . The crosses are the 
numerical solution, and the solid lines are the empirical 
closed-form expression in Eq. (48). 


the interface is infinitely stiff and the static friction law 
approaches AC’s friction law. As k decreases the front 
velocity also decreases and v c —>• I in the k —>• 0 limit 
where the interface is infinitely soft. 

The rupture criterion, Eq. (43), is essentially a crite¬ 
rion for the displacement of a block relative to its at¬ 
tachment point. For a given pre-stress r, as the interface 
stiffness k is reduced, blocks move a larger distance before 
detaching. In the limit k —> 0 the rupture front will es¬ 
sentially become a displacement wave which moves with 
a velocity equal to the velocity of sound. This explains 
the behaviour of the model seen in the k —> 0 limit. 

Unfortunately we have not been able to obtain an an¬ 
alytical solution for the front velocity in the model with 
a tangential stiffness of the interface. Instead we have 
found the empirical expression 

f(v c ,k) = ( 1 —n-" 1 ) 17 " 2 , (48) 

where the coefficients n\ and n-i are functions of k. to 
yield satisfactory agreement with the numerical solu¬ 
tions. Best fit values of the coefficients n\ and 712, for 
the values of k in Fig. 5, are given in Table I. These 
were obtained by solving the steady-state equations (Ap¬ 
pendix A 2) numerically as described in Appendix B 2 
and fitting with Eq. (48) using a least squares method. 
The predictions made by Eq. (48) are shown as solid lines 
in Fig. 5. 

C. Behaviour of the complete model 

Figure 6 shows the evolution of front velocity as a func¬ 
tion of pre-stress in the complete model, i.e. where both a 
relative viscosity and a tangential stiffness of the interface 
are included. Several values of fj and k are used, demon- 


k 

m 

n 2 

10 3 

2 

2 

10 1 

3.75 

1.7 

10° 

7.86 

1.97 

nr 1 

20.8 

2.01 

to -3 

106 

2.87 


TABLE I: Best fit coefficients to be used in the 
empirical expression for the front velocity given in 
Eq. (48). Note that the empirical approximation for 
k = 10 3 is identical to Eq. (30). 



FIG. 6: Colour online. Steady-state front velocity v c as 
a function of initial shear stress, f, for several different 
values of relative viscosity fj and interface stiffness k. 


strating that the complete model behaves in a way quali¬ 
tatively consistent with the results of Section III A 2 and 
Section IIIB. In particular, front speed increases both 
with increasing fj and with increasing k. 


IV. DISCUSSION 

We have found that, in general, the rupture front veloc¬ 
ity increases with increasing pre-stress in our ID spring- 
block-like models of extended frictional interfaces, as seen 
in Figs. 3 to 6. This is in agreement with observations 
on poly(methyl methacrylate) interfaces by Ben-David 
et al. [5]. We have also found the governing pre-stress 
parameter to be 


Ms - Mk 

That is, in addition to the pre-stress r itself, the steady- 
state rupture velocity also depends on the local friction 
parameters p s and and on the applied normal load 
through p. The same parameter f has successfully been 
applied to 2D spring-block models to scale non-steady- 
state front velocities obtained with different model pa¬ 
rameters (see Fig. 4b in [9], or [13]). It is also equivalent 
to the S ratio used in the geophysical literature [38], de- 
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fined as 


°V ~ 

(J\- <J{ 


(50) 


where <r y , <J\ and oy are the yield stress, initial stress 
and sliding frictional stress, respectively. In terms of the 
parameters in the model discussed here, a y = /i s p, ay = 
r, erf = /KkP, and we can express S' in terms of f: 


Ms - r/p __ 1 

r/p - Mk r 


(51) 


The parameter f is therefore much more general than the 
derivation from the present model alone would indicate. 

The second parameter of importance to steady-state 
ruptures is the relative viscosity parameter fj discussed 
in Section III A 2. It provides a simple way of introduc¬ 
ing energy dissipation that will occur during deformation 
of the slider and also removes resolution-dependent oscil¬ 
lations in Burridge-Knopoff-like models [12, 39], with the 
recommended value fj = y/OA. Note that the viscosity rj 
considered here is a bulk viscosity affecting the relative 
motion of blocks. It is thus qualitatively different from 
an interfacial viscosity that would affect the absolute mo¬ 
tion of a single block on the track, as was sometimes 
introduced at the micro-junction level in multi-junction 
models (see e.g. [21]) or directly at mesoscales as a ve¬ 
locity strengthening branch of the friction law at large 
slip velocities (see e.g. [40]). In our case, the friction 
force on a block in the sliding state is /lyjp. independent 
of slip speed. At a given f, increasing fj increases the 
steady-state front velocity, see Fig. 4. As discussed in 
Section III A 2 this is due to the added shear force aris¬ 
ing from the damping of relative motion between blocks. 
The particular choice fj = VOA (used in [9, 12]) is in 
Fig. 4 seen to only modestly increase the front velocity 
compared to fj = 0. 

The third and last parameter studied here is k = Ay/fc, 
the interface to bulk stiffness ratio, discussed in Sec¬ 
tion HIES. In Fig. 5 the limit k —> oo is seen to yield 
Amontons-Coulomb-like behaviour, while decreasing k 
yields decreasing front velocities. In fact, as k —► 0 the 
front velocity approaches the speed of sound as discussed 
earlier due to the front becoming a sound wave. These 
results should be relevant to various similar ID and 2D 
models in which blocks are elastically connected to the 
base by springs [9, 11-13, 21, 22]. 

Recent simulations, in a 2D spring-block model with a 
friction law at the block scale emerging as the collective 
behaviour of many micro-junctions in parallel, have iden¬ 
tified two different slip regimes for individual blocks [11]. 
A fast (inertial) slip regime is followed by a slow slip 
regime, controlled by the healing dynamics of the inter¬ 
face after rupture. Fronts driven by fast slip are fast 
inertial fronts, whereas fronts propagating when a signif¬ 
icant portion of the slipping blocks are in the slow regime 
are slow [13]. In this context, all fronts observed in the 
present ID models are of the fast type. 


Although as seen in Fig. 2a the transient front velocity 
is often sub-sonic, in our model, all steady-state fronts 
are supersonic, i.e. v c > 1. The fronts can propagate at 
arbitrarily large speeds as long as the pre-stress f is large 
enough. This has been discussed previously by Knopoff 
[34]. The velocity of sound in a ID model is the longitu¬ 
dinal wave speed, while shear and Rayleigh waves do not 
exist. Nevertheless, we think it useful to point out that 
super-shear fronts have recently been observed in model 
experiments [5], and that in the geophysical community, 
such fronts have been both predicted theoretically [41, 42] 
and confirmed experimentally [1, 6, 43, 44]. 

As seen in Fig. 2b, and previously found in [13], the ini¬ 
tial transient in front speed before steady state is reached 
can be very long when f is close to zero. Because the di¬ 
mensionless equations of motion do not change with the 
model resolution, it is clear that in this model, the length 
of the transients is given by a fixed number of blocks 
rather than a physical length, so we have overcome the 
problem of getting close to the steady state by perform¬ 
ing simulations with a large number of blocks. In the 
other limit of f —► 1, the transient length vanishes. We 
provide a demonstration of this result in Appendix C. 

To investigate the transient length’s dependence on f 
for small values of f we initialise the system with a con¬ 
stant pre-stress as in Fig. 3 and let the rupture propa¬ 
gate until its velocity has reached 97 % of the steady-state 
value. We define the point at which this happens as the 
transient length and plot in Fig. 2b the transient length 
as a function of f. For small values of f the transient 
length is very large. 

The above considerations emphasise the fact that in 
spring-block models like the one studied here it is impor¬ 
tant to choose the resolution carefully: this choice will 
indeed select the physical length of transients in front 
dynamics. The size of each block should be equal to the 
screening length [45], which in a purely elastic model is 
given by A d ~ d 2 /a , where d is the distance between 
micro-contacts and a is the size of micro-contacts. For 
micrometer-ranged roughnesses, we expect a ~ 1 pm and 
d ~ 10 pm to 100 pm, yielding A ~ 0.1 mm to 10 mm, i.e. 
in the millimeter range. The typical horizontal length 
scale for extended lab-scale interfaces is L ~ 100 mm, 
which yields N ~ 100, which is consistent with the num¬ 
ber of blocks used here or in our previous studies [9,11 

13]- 

Let us now compare our results to those of previ¬ 
ous studies of steady state rupture velocities. In the 
Amontons-Coulomb case, our solution v c (t) takes the 
exact same form as the one found in [32] for the Burridge- 
Knopoff model with Amontons-Coulomb friction (com¬ 
pare Eq. (B8) and Eq. (A14) in [32]). As in [32], we find a 
well-defined continuum limit where rupture velocities do 
not depend on the chosen resolution. Compared to the 
studies in [30, 31] of the Burridge-Knopoff model with 
velocity weakening friction, our results are qualitatively, 
although not quantitatively, similar. In particular, we 
also find that the rupture velocity increases with increas- 
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ing shear prestress of the interface and with increasing 
values of the viscous coefficient. Note that [30-32] did 
not discuss the effect of an interfacial stiffness on rup¬ 
ture speed. 

Trpmborg et al. [11], Trpmborg et al. [13] showed that, 
in their 2D model, the rupture and slip velocities are pro¬ 
portional. In ID it is possible to derive a similar, exact 
relationship between the average slip speed and rupture 
velocity in the case of Amontons-Coulomb friction and 
no viscosity. This derivation is provided in Appendix D, 
where it is found that the local rupture velocity v Ct i at 
block i in a simulation is related to the local average slip 
velocity u Ijavg by 


v 


c,i 


avg 

1 — f 


(52) 


This relationship is demonstrated in Fig. 8 , where we 
have plotted the local rupture velocity as a function 
of the local average slip velocity for various values of 
f as measured in simulations similar to those seen in 
Fig. 2a. Rescaling the average slip velocity by 1/(1 — f), 
a straight line of unit slope is obtained. For com¬ 
parison with the rescaling formula found by Trpmborg 
et al. [11], Trpmborg et al. [13], it is instructive to write 
Eq. (52) using dimensional quantities. Equations (10), 
(18) and ( 22 ) yield 

_ <avg£S 
Vc.i ? 

fl s p -T /i s P-T 

where we have used a = L/{N— 1) and k = (N— 1 )ES/L 
as in [12], where L is the length of the slider, E is Young’s 
modulus and S is the cross-sectional area of the slider. 
This strongly resembles the rescaling formula in [11, 13] 
which has the same denominator, while the characteristic 
force in the numerator is different due to the different 
interfacial laws applied. 

An important question is whether the results obtained 
in the present ID model can be extended to 2D models. 
To answer this question, we perform a series of simula¬ 
tions using the 2D model described in [11, 13], with model 
parameters suitable for the study of steady state front 
propagation. In particular, the slider’s length is twenty 
times larger than in [ 11 ], so that front propagation has a 
chance to converge towards a steady state (the length of 
transients ranges from less than 40 blocks for f = 0.95 
to longer than the system length for f = 0.2). We first 
choose a reference set of parameters, in which parame¬ 
ters are the same as in Table SI of [11], except L = 2.8 m, 
N x = 1140, M = 1.5 kg, F n = 3840 N, f tr i gge r = 0.95, 
and a = 77 /I 6 O. The rest of the settings are as follows: 
The width of the initial junction force distribution is zero, 
so the interface springs effectively act as a single spring 
per block. We checked that the mean junction slipping 
time tR, while important for slow fronts, does not affect 
these fast front results. Steady-state front velocity is es¬ 
timated as in Section III A 1, by fitting a straight line to 
v c (l/n) for the last (~ N x / 2) blocks towards the leading 
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FIG. 7: Colour online. Steady-state front velocity, v c . 
as a function of initial shear stress, r, for simulations 
using the full 2D model of [11], for several different 
values of mass, damping coefficient, bulk stiffness and 
interface stiffness. Five sets of parameters were used: 
(i) Reference set; (ii) mass doubled and damping 
coefficient increased by a factor of \/2, leaving fj 
unchanged; (iii) bulk and interface stiffnesses halved, 
leaving k unchanged, (iv) damping coefficient and fj 
reduced by a factor of 4; and (v) interface stiffness and 
k increased by a factor of 10. The dashed line plots 
Eq. (36) for the reference value of fj. Front speeds are 
normalized against the longitudinal bulk wave speed. 
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edge (excluding the last ~ 40, which have v c increasing 
due to an edge effect). 

We first run simulations of the reference model for vari¬ 
ous values of r and plot the normalized steady state front 
speed v c as a function of f in Fig. 7 (blue crosses). We 
observe that the results are qualitatively fully consistent 
with the behaviour in ID, i.e. Fig. 6. In particular, 
the steady state front speed is always supersonic, tends 
towards the longitudinal bulk wave speed for small pre¬ 
stress and diverges for large prestress. 

We then show that the two control parameters iden¬ 
tified in ID, fj and k, are also controlling the 2D front 
speed. To do this, we change model parameters (slider’s 
mass, damping coefficient, bulk and interfacial stiff¬ 
nesses) in such a way that the rescaled parameters fj and 
k are kept unchanged. Figure 7 clearly shows (see cyan 
dots and red squares) that these changes do not affect 
the values of v c . 

Finally, we show that changes to the values of either fj 
or k induce variations in the v c (f) curve which are fully 
consistent with those observed in ID (Fig. 6): decreasing 
fj decreases the front speed (see black diamonds) whereas 
increasing k increases the front speed (green stars). All 
these results indicate that our main findings from the ID 
model are not specific to ID, but also hold in 2D, with 
the same non-dimensional control parameters and all the 
same qualitative characteristics. 
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V. CONCLUSIONS 

We have systematically studied the quantitative de¬ 
pendence of the steady-state rupture front velocity on 
pre-stress, damping and stiffness of the interface in a ID 
spring-block model. We find that front velocity changes 
significantly by changing any of these parameters, the 
result of which can be seen in Figs. 3 to 6. 

Increasing pre-stress leads to increasing rupture veloc¬ 
ities, in agreement with both experiments [5] and 2D 
models (see discussion and [9, 10, 13]). Specifically, 
for the model with no viscosity and no finite stiffness of 
the interface, we derive a closed-form expression for the 
front velocity, given by Eq. (30). The dimensionless pre¬ 
stress parameter found to be controlling front velocities, 
f = {r/p — Hk)/(n s — p, k), also depends on the strength 
of the interface through the frictional parameters /x s and 
/Jk. It is essentially a version of the S ratio often used by 
the geophysical community [38] and shown to also apply 
in 2D models [9, 13]. 

Material dampin g affe cts the front velocity through the 
parameter fj = p/y/km. Increasing values of fj are seen 
to yield increasing front velocities caused by the addi¬ 
tional shear force. A semi-empirical expression for the 
dependence of the front velocity on fj and f is given in 
Eq. (36). 

Front velocities are seen to decrease with decreasing 
tangential stiffness of the interface through the parame¬ 
ter k = kt/k, i.e. the ratio between the shear stiffness 
of the interface and the material stiffness of the slider. 
An empirical expression for this dependence is given in 
Eq. (48). In fact, in the limit of a very soft interface 
compared to the material stiffness, steady-state rupture 
velocities are seen to approach the velocity of sound. The 
qualitative behaviour of all these parameters are seen to 
carry over to a model with both a finite stiffness of the 
interface and relative viscosity. 

From Fig. 2b it is clear that transients can become 
very long, especially for low pre-stresses. This, coupled 
with a heterogeneous interface where t can be negative, 
suggests that experimentally observed rupture fronts like 
those in e.g. [4, 5] may be dominated by transients. Di¬ 
rect comparison between these rupture fronts and the 
ones studied here may therefore not be possible, but the 
qualitative behaviour on parameters such as the internal 
damping and interface stiffness should remain valid. 

Also note that viscoelastic materials have been shown 
in finite element simulations to exhibit memory effects. 
Stress concentrations left at the arrest location of one 
precursory slip event are not erased by the following rup¬ 
ture [46]. This causes non-homogeneous initial stresses 
for subsequent events. An interesting direction for fu¬ 
ture work would be to investigate the transient speeds 
resulting from such complicated stress states. 

Despite the limitations of the model discussed above, 
it can provide valuable insight into rupture dynamics of 
frictional interfaces. We have shown here that the ID 
results can be extended to a 2D model with the same in¬ 


terfacial law and bulk damping. Experimentally, it would 
be interesting to investigate the dependence of front ve¬ 
locities on the interface stiffness and bulk viscosity as 
they have been shown here to affect the rupture velocity 
significantly. Also, investigating the dependence of the 
front velocity on the system length would shed light on 
the influence of transients on observed ruptures. 
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Appendix A: Deriving the equations for the 
steady-state rupture velocity 

Here we derive the equations for steady-state rupture 
in the models discussed in this paper. We assume u n (t) > 
0 for all detached blocks, in agreement with the numerical 
solution of the equations for all cases where f + > 0, as 
discussed in Section III. The system is considered to be 
infinitely long and the tangential driving velocity much 
smaller than the front velocity, we can therefore ignore 
the system boundaries. 

1. Amontons—Coulomb friction 

Here we derive the equations for the steady-state front 
velocity for the model with AC friction. As our start¬ 
ing point we use the dimensionless equations of motion, 
Eqs. (11) to (13), (21) and (25). Consequently, the con¬ 
trolling parameters in the equations of motion are f = f + 
and fj. 

The equations of motion for moving blocks are given 

by 

Un = ^n-f-1 2 U n T U n —\ (Al) 

+ fj(u n+ i - 2u n + M n _i) + f. 

To eliminate the dependence on f we introduce u n de¬ 
fined by 

u n =f ( u n + P/2) , (A2) 

where the acceleration of blocks due to the force f is taken 
into account explicitly by the term fP/2 in Eq. (A2). 
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Equation (Al) simplifies to 

Vn — &n +1 2'U n T U n — i 

_ . . (Aoj 

+ fj(u n +i - 2 u n + U n - 1), 

which is our final equation of motion for moving blocks. 
If the front is propagating at a constant velocity, 


which relates the solution of Eqs. (A3), (A5), (A6) 
and (A8) for a given v c to the corresponding f. Note 
that these equations, in the case of fj = 0, reduce to the 
equations determining the front velocity in the Burridge- 
Knopoff model using AC friction and the approximations 
of slow and soft tangential loading as derived by Muratov 

[32]. 


u n (t) = u n+1 (t + At), u n (t) = u n+ i(t + At), (A4) 

must hold, where At is the time between the triggering 
of two neighbouring blocks. It is therefore sufficient to 
consider the system in a time interval of length At. We 
choose t € [ ti,ti + At], where block i begins to move at 
t = ti- For convenience, and without loss of generality, 
we choose t, = 0. Using Eq. (A2) this yields 

fii(0) = 0, fii(0) = 0, (A5) 


2. Including a finite tangential stiffness of the 
interface 

Here we derive the steady-state equations for the model 
including a tangential stiffness of the interface discussed 
in Section IIIB. The equations of motion for sliding 
blocks are given by Eq. (A3) as their equation of mo¬ 
tion is identical to the AC case. For stuck blocks we 
combine Eqs. (11) to (13), (41) and (A2) and get 


which is the initial condition for block i. For convenience 
we choose i = 0. 

The initial condition for Eq. (A3) can be obtained by 
evaluating Eq. (A4) at t = 0 and using Eq. (A2). This 
yields 


«n(°) = u n +i (v c *) + A, 

1 (A6) 
&n(0) = U n+ 1 (v c 1 ) + — , 

V c 

since 

A t = = u- 1 (A7) 

V c 

from the definition of X, Eq. (10). 

The equation of motion for all moving blocks is given 
by Eq. (A3), but the equation of motion for block i 
can be rewritten taking into account that block i + 1 
is stationary. Inserting u,; + i = 0 into Eq. (A2) yields 
u i+1 = —P/2, and using Eq. (A3) with n = i we have 


Ui = Ui-i — 2ui — P/2 

+ fj{ui -1 - 2 Hi - t), 


(A8) 


Vn — Wt +1 T Vj n — i 2u n 

+ fj(u n+1 - 2u n + u n -i) (All) 

— k(u n + P/2) — 1, 

where, without loss of generality, we have set u® tlck = 0. 
For convenience we again let block * = 0 detach at t = 0, 
and by the same argument as above these equations need 
only to be solved for t £ [0, At]. The initial conditions are 
the same as for the AC case, Eq. (A6), with the exception 
that only blocks far away from the rupture front keep a 
constant position. The initial conditions are therefore 
given by 


Vn— 1(0) — V n [v c ) T 

2 v c 

U„_l(0) = UniVc 1 ) + —, 

V c 


(A12) 


'U j n—> oo — 'U , n—> oo — 0* 

Using Eqs. (43) and (A2) at t = 0 with u^ lck = 0 we get 
ku n ( 0) = _ T , (A13) 


Solving Eqs. (A3), (A5), (A6) and (A8) result in u n (t) 
for a given rupture velocity v c . This velocity is related to 
the parameter t through the rupture criterion, Eq. (24), 
at time t = 0 for block i = 0. At t = 0 we have ito(0) = 
ui(0) = uq (0) = ui(0) = 0, and using Eqs. (12), (13), 
(24) and (A2) in addition to t = /+ we get 

V— i(0) +fju- 1(0) = 1 _ T , (A9) 

T 

or equivalently 


U-l(0) + TJU- 1(0) + 1 


which relates the solution of Eqs. (A3), (All) and (A12) 
for a given v c to f. 


Appendix B: Solving the equations for the 
steady-state rupture velocity 

Here we solve for the steady-state velocity analytically 
using the equation set derived in Appendix A1 for the 
model with AC friction. In addition we describe a numer¬ 
ical solution procedure that we use to solve the steady- 
state equations for the models with either AC or elasto- 
plastic friction laws. 
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1. Analytical solution of the Amontons—Coulomb 
steady-state rupture velocity equations 

We solve Eqs. (A3), (A5), (A6) and (A8) using the 
iterative approach employed by Muratov [32] where the 
solution is a power expansion in 



If the front velocity is large, i.e. in the limit z = v~ l —» 0, 
the distance required for block 0 to move in order to ini¬ 
tiate movement of block 1 is negligible. Therefore, block 
0 is stationary in this limit. Also, when the rupture ve¬ 
locity is infinitely high, interactions between the blocks 
become less prominent because the time interval At be¬ 
comes negligibly small. 

To zeroth order in z, we therefore ignore the interaction 
terms in Eqs. (A3) and (A8), which yields «n°' ) = 0, i.e. 

«£» =<*£»*+ &<?>, (B2) 

where a^°' ) and are constants to be determined from 
the initial condition. Using Eq. (A6) we get the two 


coupled difference equations 

a n ] = 4+1 + U ( B3 ) 

4 0) = 4+i* + 4+i + ( B4 ) 

where = b 4 = 0 from Eq. (A5). The solution is 


4 0) = ~ nz , 

(B5) 

6 (0) _ 4.2 

On 2 Z ’ 

(B6) 


and we therefore have 

n 2 

= — nzt + — z 2 , n A 0. (B7) 

This is the first iteration, giving the zeroth order solu¬ 
tion to Eqs. (A3), (A5), (A6) and (A8). The zeroth or¬ 
der relationship between v c and f is obtained by using 
Eqs. (A10) and (B7). 

The first order solution is obtained by substituting u n 
and u n using Eq. (B7) into Eqs. (A3) and (A8), inte¬ 
grating twice with respect to i and then using Eqs. (A5) 
and (A6) to eliminate the integration constants. This 
approach rapidly becomes cumbersome, so we have used 
Mathematica to go to higher orders. Here we simply 
give the solution: 


z 2 

— (r?-l)(r? +1) + 

2 3 _ _ 

gj-»?(»7- 1 )(V+ 1) - 

- l){fj + l)(7fj 2 - 3) — 

^v(v - 1 )(V + 1 )( 19 ? 7 2 - 11 ) + 

.6 

— (fj - 1 )(fj + 1)(229?7 4 - 226 fj 2 + 45) + 

6! 

z 7 

—fj(fj - 1 )(fj + 1)(995?7 4 - 1154f7 2 + 295) - 

.8 

— (fj - 1 )(fj+ 1)(1715177 6 - 26837^ 4 + 12349?7 2 - 1575) - 
8! 

z 9 

—fj(fj - l)(fj + l)(102083r? 6 - 178177?7 4 + 95033^ 2 - 14971) + 

9! 

z 10 

— (fj - l)(fj + 1)(229414177 s - 49303847? 6 + 3640514?7 4 - 10638167 2 + 99225) + 
0(z 41 ). 


(B8) 


We conclude this section with a brief discussion of the validity of the solution in Eq. (B8). It is valid only for 
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f < 1 due to the form of the rupture criterion. It is not 
valid for f = 0 because of the transformation in Eq. (A2). 
If f < 0, then u n becomes negative (combine Eqs. (A2) 
and (B7)). As a result, the velocity must be negative, 
but it was assumed to be positive. Equation (B8) is con¬ 
sequently not valid for t < 0, i.e., we get the constraints 

0 < t < 1, (B9) 

which is consistent with the numerical results in Sec¬ 
tion III. 


2. Numerical solution procedure 

To solve the steady-state equations derived in Ap¬ 
pendix A numerically, we used the following iterative 
scheme: 

1. Select the desired z = v~ x for which the value of 
f is to be found. An initial guess for the position 
and velocity of all blocks must be made at t = 0, 
we used u„(0) = u n (0) = 0 as the initial guess for 
all results presented here. 

2. Solve Equations (A3) and (A8) or Eqs. (A3) 
and (All) using a numerical scheme for solving 
differential equations, e.g. the fourth order Runge- 
Kutta method. 

3. Calculate a new estimate for the initial conditions 
using Eqs. (A5) and (A6) or Eq. (A12) and the cho¬ 
sen front velocity. In the model with a tangential 
stiffness of the interface we have u n ^oo = Mn->oo = 
0, which is implemented as un = Un = 0 where N 
is the number of blocks in the calculation. We find 
that N = 100 yields satisfactory results. Calculate 
f using Eq. (A10) or Eq. (A13). 

4. Repeat steps 2. and 3. until t has converged. The 
solution has converged when the difference in t be¬ 
tween two iterations is less than some tolerance e. 
We use e = 1CT 6 . 

The functions v c (f) or equivalently f(v c ) can be calcu¬ 
lated by repeating the above steps for several values of 

v c . 


Appendix C: Deriving the r — > 1 limit of transient 
length 

Here we show that in the limit of f —> 1 the transient 
length vanishes by considering the motion of the trailing 
edge block as it begins to move. 

Using Eqs. (11) to (15) with AC friction [Eq. (26)] the 
equation of motion for the trailing edge block becomes 

ute = — wte + t + .Ft (Cl) 

where we have set ?7 = 0 for simplicity. As before we 
assume that FY is independent of time. The initial con¬ 
dition for the trailing edge block is wte(£ = 0) = ute(£" = 


0) = 0 and rupture initiates when t + FY = 1. This yields 

Ute© = 1 — cos t. (C2) 

The rupture criterion, Eq. (24), applied to the second 
block from the trailing edge yields 


wte(U) = 1 — cosfi = 1 — r, (C3) 

where t\ is the time at which the motion of the second 
block from the trailing edge is triggered. Noting that 
z\ = ti = l/v c ,i, we have 

2 4 

T = cosz = 1 - + A- + 0(z 6 ). (C4) 

Comparing with Eq. (28), the series of the steady-state 
rupture velocity, the first two terms are identical, i.e. 
as f —»• 1 (z —> 0), the rupture will instantly reach the 
steady-state velocity and the transient length approaches 
0. For smaller values of f (larger values of z), the tran¬ 
sient length increases due to the deviations in the higher 
order terms. 


Appendix D: Slip speed vs rupture speed 


Here we derive the relationship between the rupture 
and slip velocity in the model with Amontons-Coulomb 
friction and no viscosity. Consider block i, which started 
to move at time t = U, causing an increased shear force 
on block i + 1. We only consider blocks moving in the 
positive direction, i.e. only ff = f is required. The local 
rupture velocity is given by Eq. (A7), 

^ c ’ i= Al? (D1) 

where U = U+i — U is the time between the triggering of 
block i and i + 1. The position of block i at t, is Ui = 0. 
The increase in shear force required for block i + 1 to 
start moving is given by the rupture criterion, Eq. (24), 
and the position of block i at U+i is therefore given by 


Ui{t i+ 1) = 1 - T. 

The average slip velocity is consequently 

^ _ Ui(ti+ 1) - Ui{ti) _ 1 - r 

Uj ’ avs ^ A U ~ A U 


(D2) 


(D3) 


and the rupture and slip velocities are related by 


v c ,i = 


'U'i, avg 

1 — f 


(52) 


This relationship is demonstrated in Fig. 8, where we 
have plotted the local rupture velocity as a function of the 
local average slip velocity for various values of t as mea¬ 
sured in simulations (similar to the simulations shown in 
Fig. 2a). Rescaling the average slip velocity by 1/(1 — f), 
a straight line with unit slope is obtained. 
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(a) Local rupture velocity v c as a function of local slip velocity 

^avg- 



(b) Local rupture velocity v c as a function of the rescaled local 
slip velocity u aV g/(l — f). The black line is a straight line through 
the origin with a slope of unity as predicted by Eq. (52). 



FIG. 8: Colour online. Rupture velocity and slip velocity are closely related. These results have been obtained using 

simulations as in Fig. 2a for f = 0.1, 0.2,..., 0.9. 
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